Predicting the critical density of topological defects in 0(N) scalar field theories. 
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O(N) symmetric Xcj> 4 field theories describe many critical phenomena in the laboratory and in the 
early Universe. Given N and D < 3, the dimension of space, these models exhibit topological defect 
classical solutions that in some cases fully determine their critical behavior. For N — 2, D = 3 it 
has been observed that the defect density is seemingly a universal quantity at T c . We prove this 
conjecture and show how to predict its value based on the universal critical exponents of the field 
theory. Analogously, for general TV and D we predict the universal critical densities of domain walls 
and monopoles, for which no detailed thermodynamic study exists. This procedure can also be 
inverted, producing an algorithm for generating typical defect networks at criticality, in contrast to 
the canonical procedure 111, which applies only in the unphysical limit of infinite temperature. 
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0(N) symmetric scalar field theories are a class of 
models describing the critical behavior of an great variety 
of important physical systems. For example, for N = 3 
they describe ferromagnets, the liquid vapor transition 
and binary mixtures for N — 1 and superfluid 4 He and 
the statistical properties of polymers, for N = 2. In the 
early Universe N — 2 describes the phase transition as- 
sociated with the breakdown of Peccei-Quinn symmetry, 
and models of high energy particle physics may belong 
to the universality class of O(N) scalar models, whenever 
the mass of the Higgs bosons is larger that that of the 
gauge bosons. O(N) scalar models are also invoked in 
most implementations of cosmological inflation. 

One of the fundamental properties of O(N) X\<j>\ 4 field 
theories is the existence, for N < D, of static non-linear 
classical solutions, (domain walls, vortices, monopoles) 
that we will refer to henceforth as topological defects. 
At sufficiently high temperatures, topological defects can 
be excited as non-perturbative fluctuations. Their dom- 
inance over the thermodynamics, due to their large con- 
figurational entropy, is known to trigger the phase tran- 
sition in 0(2) in 3D and 2D, and their persistence at low 
energies prevents the onset of long range order in 0(2) 
D < 2 and in 0(1) in ID. 

It is therefore natural that the universal critical expo- 
nents characterizing the phase transition in terms of de- 
fects and through the behavior of field correlators must 
be connected. This connection is made more quantitative 
whenever one can construct dual models, field theories 
which possess these collective solutions as their funda- 
mental excitations || . In the absence of supersymmetry 
rigorous mappings between the fundamental models and 
their dual counterparts exist only in very special cases 
|^],[|. Duality has been suggested and empirically ob- 
served to be a much more general phenomenon, though. 

In this letter we explore the duality between the criti- 
cal behavior of the 2-point field correlation function and 



defect densities at criticality. We will show that it leads 
to the proof that the critical density of vortex strings, ob- 
served in recent non-perturbative thermodynamic studies 
of 0(2), is a universal number. Among other insights Q] 
this shows that the phase transition in 0(2) in 3D occurs 
when a critical density of defects is reached, connecting 
directly the familiar picture of the Hagedorn transition 
in vortex densities to the more abstract critical behavior 
of the fields. We also extend our procedure to different N 
and D, making predictions for the values of the universal 
densities of domain walls and monopoles, in 2 and 3D. 

Finally the inversion of this procedure allows us to eas- 
ily generate typical field configurations at criticality. This 
is of fundamental practical importance. Recent experi- 
ments in 3 He and large scale numerical studies of the 
theory |(| have lent quantitative support to the ideas, due 
to Kibble and Zurek |8| , that defects form at a second 
order phase transition due to critical slowing down of the 
fields response over large length scales, in the vicinity of 
the critical point. Defect networks hence formed have 
densities and length distributions set by thermal equilib- 
rium at T = T c + . 

In contrast most realizations of defect networks used, 
eg., in cosmological studies are generated using the 
Vachaspati-Vilenkin Q] (VV) algorithm. It relies on lay- 
ing down random field phases on a lattice and searching 
for their integer windings along closed paths. The abso- 
lute randomness of the phases corresponds to the T — > oo 
limit of the theory. More fundamentally it yields defect 
networks that are quantitatively distinct from those in 
equilibrium at criticality, i.e., at formation. 

Fig. |l| shows the behavior of a system of vortex strings 
at a second-order phase-transition, for 0(2) in 3D. The 
data was obtained from the study of the non-perturbative 
thermodynamics of the field theory At T c the total 
density of string p t ot displays a discontinuity in its deriva- 
tive, signaling a second order phase transition. 
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A disorder parameter can be constructed in terms of 
string quantities by dividing the string population into 
long string (typically string longer than ~ L 2 , where L 
is the size of the computational domain) and loops, com- 
prising of shorter strings. The corresponding densities 
are denoted by piong and pi oop . In Fig. Qwe can observe 
that piong consistently vanishes below T c , except for a 
small range of /3 where it increases rapidly to a finite 
critical value. In || we conjectured that in the infinite 
volume limit p; n f exhibits a discontinuous transition. 
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FIG. 1. The string densities, total, loops and long string, as 
a function of inverse temperature /3. At j3 c , the densities dis- 
play derivative discontinuities, signaling a second order phase 
transition. ptot (/3c) coincides for different studies, leading to 
the conjecture that it is a universal number. 

The value of the total string density at /3 C , Ptot(/?c) — 
0.20 coincides with results from studies of different mod- 
els in the same universality class ||,[l0| . This fact lead us 
to the conjecture Q that ptot(Pc) is universal. 

In order to prove this conjecture we appeal to a 
well known result, due to Halperin and Mazenko [ pd| . 
Halperin's formula expresses po, the density of zeros of 
a Gaussian field distribution in terms of its two-point 
function. For an O(N) theory the relevant quantity 
is the O(N) symmetric correlation function G(x) = 
(</>(O)0(a;)' 1 '), resulting in 
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Eq. (|l]) measures the density of coincident zeros of all 
iV components of the field at a point. Coincident zeros 
occur at the core of topological defects. Depending on 
N and D, coincident zeros can be interpreted as either 
monopoles, strings or domain walls. In the particular 
case of a Gaussian 0(2) theory in D — 3, Halperin's 



formula allows us to compute the density of vortex strings 
crossing an arbitrary plane in three dimensional space, a 
quantity that is clearly proportional to p to t- 

The last key observation is that in the critical domain 
of a second order transition, all O(N) theories are ef- 
fectively approximately Gaussian, but with non-trivial 
critical exponents. In particular renormalization group 
analysis shows that the mass and quartic coupling vanish 
at T c Higher order polynomial terms (eg. oc 6 ) 

may be generated but are small. Hence in the critical 
domain the field two-point function can be written as 
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where r\ << 1. is the universal critical exponent taking 
into account deviations from the mean- field result. 

Thus the effective Gaussianity of the theory allows us 
to use Halperin's result to compute the critical value of 
Ptot (/?c)- Note that, modulo renormalization, the final re- 
sult depends only on r\ establishing, as conjectured, that 
Ptot(/?c) is a universal quantity. 

Substituting, Eq. (§) into @) we obtain: 
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where we have introduced upper and lower momentum 
cut-offs fc max and fc m ; n . In the case of a lattice of size L 
and lattice spacing a we take k mul = 2n/L and fc max = 
27r/a which leads to 
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For large enough lattices, a/L « 1, and we obtain 
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In order to generate quantitative predictions we need 
to determine the exact proportionality factor in Eq.(^|). 
This can be achieved by invoking the other instance when 
the interacting theory becomes Gaussian. In the high 
temperature limit (3 — > 0, the effective interaction be- 
comes irrelevant. In terms of renormalization group the 
theory displays a (trivial) Gaussian fixed point, with van- 
ishing effective coupling. On the lattice fields at different 
points will be completely uncorrelated. Note that this 
situation corresponds to the usual VV algorithm where 
a field is thrown randomly on the lattice and a network 
of strings is built by identifying phase windings. Fig. |l| 
shows the agreement of the densities observed at = 0. 
with the well known VV result of ptot = 1/3 with 75% 
long string. Since the totally uncorrelated field corre- 
sponds to a flat power spectrum G(k) ~ k° we normalize 
Halperin's expression by imposing ptot = 1/3 for r) = 2, 
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rj is always much smaller than 1. Setting r\ = we obtain 
p tot {T c ) = 5/27 = 0.185, close to the value p tot ^ 0.2 
observed both in the A</> 4 §] and XY @ studies in 3D. 
A more accurate estimate can be obtained by replacing 
77 by its precise value for the universality class to which 
these theories belong. From Jl2| we have r\ ~ 0.035 we 
get ptot = 0.190, closer to the observed value. 

A similar exercise permits the computation of the crit- 
ical density of domain walls for O(l) and monopoles in a 
0(3) theory at the critical temperature in 3D. 

The density of domain walls per link is 
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Note the value of ptot(P = 0) = 1/2 corresponding to the 
high-temperature limit. At /3 C we get, with r\ = 0.034 
Jl2[ , ptot(Pc) — 0.38. For monopoles we will take for the 
flat-spectrum case ptotifi = 0) ~ 0.1. A better estimate 
can be obtained from a tetrahedral discretization of the 
sphere, resulting in ptot(P = 0) = 3/32 and 
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with rj = 0.038 Q leading to the critical value ptot(Pc) = 
0.040. Finally for domain walls in 2D, the density per link 
at p c is (using p tot (/3 = 0) = 1/2): 
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Taking rj = 0.26 @ we obtain p tot (/3 c ) = 0.24. 

The present procedure can be inverted to generate a 
typical defect network at criticality. The approximate 
Gaussianity of the field theory at T c implies that the 
statistical distribution of fields, P[c/>] is given by 
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This distribution can be sampled by generating fields as 
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where R(k) is a random number extracted from a Gaus- 
sian distribution, with zero mean and unit variance. The 
field can then be Fourier transformed to coordinate space, 
its phases identified at each site, and vortices found in the 
standard way. Since we will be willing to compare results 
from this algorithm with the ones measured in lattice 
Langevin simulations we chose to employ the exact form 
of the field correlator on the lattice: 
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We have performed several tests on the algorithm, by 
comparing it to the results of the non-perturbative ther- 
modynamics of the fields at criticality. We used lattices 



of size N? at with iVj at = 16,32,64 and 128. All results 
are averages over 50 samples obtained from independent 
random realizations. Fig. [2] shows the string densities for 
values of 77 between 0. and 0.1, including all reasonable 
values of 77 in 3D. 
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FIG. 2. The string densities from the 50 Gaussian realiza- 
tions as a function of r\ for a lattice with iVi a t = 64. Error 
bars indicate standard deviation from the mean. 

The values for the densities depend on the size of the 
lattice, converging to finite values for large Aj at . In Fig. || 
we can see the scaling of ptot with box size for two choices 
of the critical exponent, the mean field value 77 = 0. and 
the theoretical result for the 0(2) universality class in 
3D, 77 ~ 0.035. We can predict the form and the power of 
this scaling through Eq. (||). Writing a/L = l/JVi at , the 
number of points in the lattice, and expanding Eq. (|4j) in 
powers of 1/iViat we see that Halperin's result converges 
to its infinite volume limit according to: 
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To check these scalings we fitted the data of Fig. || to a 
power law of the form: 



Ptot(-Nlat) = Ptot (OO) 
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For 7/ = and 77 = 0.035 we found: 



77 = 0.0, p to t (00) = 0.1969, A = 0.3259, a = 1.060; 
77 = 0.035, p to t (00) = 0.2012, A = 0.3422, a = 1.124. 

These values of a are indeed close to 1, with a larger 
correction for 77 = 0.035 as expected from Eq. (|l3|). 

In pj for a lattice of size iVi at = 100 we measured 
Ptot (/?c) = 0.198 ± 0.004. For a Gaussian field with 77 = 
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0.035 we obtain p tot — 0.203 ± 0.003. The agreement of 
the two results is very satisfactory. 
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FIG. 3. The total string density for two values of n for 
Mat = 16, 32, 64 and 128 and respective fits to a power law. 
Statistical errors are much larger than the deviation of the 
points to the fits. 



perform a precise comparison however Halperin's formula 
should be rederived for a field theory on the lattice. De- 
spite these shortcomings Halperin's formula has the merit 
of being the only analytic way of estimating the critical 
densities of defects in theories where non-perturbative 
thermodynamic results are scarce. 

We have therefore established the connection between 
the universal critical exponent characterizing the behav- 
ior of the O(N) field 2-point correlator and the critical 
density of defects. This relation implies that defect den- 
sities at T c for a system undergoing a second order phase 
transition are universal numbers. We predicted them for 
several O(N) models in 2 and 3D. Based on these insights 
we proposed a new algorithm for generating networks of 
defects at the time of formation. In particular, we have 
shown that this algorithm reproduces accurately all the 
features of a string network in 3D at criticality. This pro- 
cedure, instead of the more usual algorithm of H should 
be used to generate typical defect networks at the time 
of their formation. 

We thank R. Durrer and W. Zurek for useful discus- 
sions. This work was partially supported by the ESF and 
by the DOE, under contract W-7405-ENG-36. 



The results for pi on g and pi oop using these two different 
methods are also in good agreement. In this case we were 
not able to find a reasonable scaling expression though, 
but this is to be expected given the arbitrariness of the 
long string definition. The results for N\ at — 100, using 
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0.035 are p\ c 



0.080 ± 0.004, p lo 



0.121 ± 



0.004. These compare well with the non-perturbative 



results pi 



ong 



0.076 ±0.005, pi 



oop 



0.120 ±0.004. Even 
more impressive is that the string length distribution at 
criticality can also be reproduced by our Gaussian field 
algorithm. This distribution can be successfully fitted to 
an expression of the form pi, 



n{l)=Ar~ l e-' 3 ' 71 . 
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The fit to the results of the Gaussian field algorithm 
shows a small variation of the parameters A, 7, and a 
for 77 ~ 0. — 0.1. a is consistently zero, reflecting the fact 
that the spectrum is always scale invariant. The value 
of 7 varies between 2.34 and 2.40. For the critical ex- 
ponent 77 = 0.035 we obtained 7 ~ 2.35. Once again 
this is in good agreement with the result from the lattice 
non-perturbative thermodynamics at T c ||, 7 ~ 2.36. 

Finally the predictions for p to t from Halperin's for- 
mula, when compared to the accuracy of the Gaussian 
algorithm seem rather poor. The expression is meant to 
apply for continuum distributions, while all other values 
of Ptot were obtained on the lattice. A straight substitu- 
tion of the lattice correlators (12) into (|l|) increases ptot 
to 0.21 from 0.19, covering our full range of results. To 
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